General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



N83-34862 


(NAS A-Tfl-85079) THB MXB 1 9 16-053/40 1 9 15-05; 

BOaST PROPEETIBS AND CONSTRAINTS ON A St 
HXNOTE BXNAEI SBCONDABY (NASA) 40 p 
UC A03/NF AO I CSCL 03A Unclas 

G3/89 41903 

IW^ 

Technical Memorandum 86079 


MXB1 91 6-053/4U1 91 5-05: 
BURST PROPERTIES AND 
CONSTRAINTS ON A 50 MINUTE 
BINARY SECONDARY 


J.H. Swank 
R.E. Taam 
N.E. White 


August 1 983 


nn 


SEP1<583 

RECEIVED 

NASASri FACIUTY 
, ACCESS DEPT. ^ 


L8 




National Aeronautics and 
Space Administration 

Goddard Space Flight Center 

Greenbelt, Maryland 20771 



MXB1916-053/4U1915..05:. BURST PROPERTIES 
AND CONSTRAINTS ON A 50 MINUTE BINARY SECONDARY 


J.H. Swank 

Laboratory for Hiflh Energy Astrophysics 
NASA/Ooddard Space Flight Center 


R*E« Taatn^ 

Department of Physics and Astronomy 
Northwestern University 


and 


N.E. Whlte^ 

Laboratory for High Energy Astrophysics 
NASA/Goddard Space Flight Center 


^Research supported In pavt Ivy NSF AST-8109826 
^ATso Dept, of Physics * Astronomy, Univ. of Maryland 


2 


ABSTRACT 

Results are presented from OSO-8 and HEAO-1 A2 observations of 34 bursts 
from the X-ray burster MXB1916-053/4U1915-05 recently discovered to show a 50 
minute binary period. While 11 bursts previously reported all had similar 
light curves, 22 observed two years later show a factor of 3 range of peak 
fluxes and decay times between 3 and 20 s. Recurrence times between 
successive bursts vary between 3 and 6 hours. We derive a ratio of steady 
flux to average burst flux of ~ 12Q. A burst observed with the HEAO-1 A2 
experiment showed an Initial temperature rise to a peak black body temperature 
of ~ 3 keV followed by the cooling typical of type I bursts. The burst was 
unusual in that the apparent projected size of a blackbody source Increased by 
a factor of 3 during the cooling phase. 

A 50 minute binary period allows many different types of mass losing 
star. We attempt to use the characteristics of the burster to discriminate 
between them and to determine If for some solution, consistent with the burst 
properties, the angular momentum loss through gravitational radiation can 
drive the mass transfer. Results of thermonuclear flash models of bursts 
Imply the secondary Is not hydrogen exhausted and that the mass transfer rate 
Is greater than gravitational radiation can drive for any of the conventional 
models that have, to date, been worked out In detail. If the peak burst 
luminosities are typical of those of other bursters (assumed to emit 
Isotropically) the Impller^ distance and the steady flux also give a high mass 
transfer. A composite model for the secondary Is marginally consistent within 
a factor of 2. Several mechanisms which could promote higher rates of mass 
transfer are discussed. 
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I. INTRODUCTION 

White and Swank (1982, hereafter Paper I) and Walter et al. (1982) 
Independently reported the discovery of absorption dips In the X-ray flux from 
4U1915-05, recurrent with a period of ~ 50 minutes. Both papers gave strong 
arguments in favor of interpreting the periodic phenomenon as the orbital 
period of the system. 4U1915-05 = MXB 1916-053 is a source of type I X-ray 
bursts, each lasting 10 s, with a recurrence time of several hours (Becker 
et al. 1977*, Lewin, Hoffman and Doty 1977). This is the first X-ray burster 
to show a confirmed X-ray periodicity. As discussed in Paper I and by Walter 
et al. (1982) there are jitters and irregularities in the dips, so that they 
are not as good a clock as partial eclipses by a companion might be expected 
to be, and a picture was developed analogous to that of cataclysmic variables 
in which a gas stream from the secondary impacts an accretion disk in a 
turbulent "hot spot", and the thickened region of the disk occults the X-ray 
source. 

The 50 minute orbital period of MXB 1916-053 is very close to the 46 
minute orbital period of the 7s X-ray pulsar 4U1 626-67, the orbit of which was 
identified from optical observations of X-ray pulsations reprocessed on the 
binary companion (Middleditch et al. 1981). Failure to find any Doppler shift 
in the X-ray pulses of 4U1 626-67 and accumulating optical identifications of 
faint counterparts to this and many other non-pulsing, non-eclipsing sources 
led Joss and Rappaport (1979), prior to the discovery by Middleditch et al. 
(1981), to propose that the secondary mass should be <0.5 M^ in 4U1626-67 and 

in the bright galactic center and globular cluster sources. is low 

enough that evolutionary mechanisms would lead in reasonable lifetimes to very 
compact binaries in which the secondary might have become degenerate (Joss, 
Avni, and Rappaport 1978). Some of these sources (Cyg X-2, Sco X-1, 



4 


4U0614+09, Aql X-1) may have periods on the order of days (Cowley, Cramp ton 
and Hutchings 1979; Cowley and Crampton 1975; Marshall and Mill It 1981; Watson 
1976). However, recent Identifications of ~ 5^ periods for 4U1 822-37 and 
4U2129-47 (Mason et al. 1980; Thorstensen et al. 1979) and ~ 4^ periods for 
4U1 636-53 and 4U1 735-44 (Pedersen, van Paradijs and Lewin 1982; McClIntock and 
Petro 1981) Imply secondaries that would populate the upper end of the 
predicted range (~ 0.5 M^). Evolution Is expected to lead to shorter periods 
for these systems and MXB1916-053 and 4U1626-67 may be examples of this. 

For orbital periods less than a few hours It has been postulated that 
gravitational radiation (GR) dominates the angular momentum losses that cause 
the orbits to decay (Faulkner 1971; Paczynski 1967). The limit on how compact 
a system will become depends on the secondary's original composition and its 
history. Orbital periods < 50 min are much shorter than the limits recently 
calculated for a secondary with standard abundances (Paczynski 1981; Paczynski 
and Slenklewicz 1981; Rappaport, Joss and Webbink 1982, hereafter RJW; Whyte 
and Eggleton 1981). Both paper I and Walter et al. (1982) discussed possible 
types of unevoTved stars that could “fit Into" a 50 minute binary with a 1.4 
Mq neutron star. A consistent picture for any of these or other possibilities 
was not addressed, although It has been pointed out (Kleboom and Verbunt 1981; 
RJW) that compact binaries with main sequence secondaries could not evolve to 
the most luminous X-ray binaries with mass transfer driven by GR alone. 

While X-ray bursts are not “standard candles" In the sense of having 
unique peak luminosities (see Lewin 1982 a,b for references and discussion) 
they still potentially offer Important Information on the nature of the 
secondary that Is not available for many bulge sources. A range of observed 
peak luminosities Is deduced from sources which can reasonably be assumed to 
cluster close to the galactic center (Inoue et al. 1980; van Paradijs 1981). 
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Although the emission mdy be modified by a scattering dominated atmosphere 
(estimates have been given by Swank, Eardley and Serlemitsos 1979; 
van Paradijs, Rybicki and Lamb 1979; van ParadIJs 1982), the spectra can be 
reasonably well described by a blackbody form and a range of effective 
projected sizes for the sources has been determined. The bursts can be 
Interpreted within the framework of the thermonuclear flash model of X-ray 
bursts to estimate the mass transfer rate and deduce the composition of the 
accreted material, and hence test the consistency of the current evolutionary 
models. 

In the next section we present new observational data of X-ray bursts 
emitted from MXB1 916-053. In Sill we deduce the parameters of the system from 
the X-ray bursts. Theoretical considerations based on binary star evolution 
will be described In SIV, and used with thermonuclear flash theory to restrict 
the parameters of the system In SV. Finally, we discuss some suggested 
avenues toward a self-consistent picture. 

II. BURST OBSERVATIONS 

The observations were made with the Goddard Space Flight Center 
proportional counter detectors on the OSO-8 and HEAO-1 satellites and the 
different sensitivities and observational constraints of the two experiments 
allow different aspects of the burst properties to be Investigated. 

MXB 1916-053 was burst active during the several days of OSO-8 observations In 
March 1976 and March 1978. Both observations provide measures of the 
recurrence times and steady flux, along with limited Indications of burst 
properties such as peak luminosities, temperatures, rise and decay times and 
burst energies. Only one short 4 hour pointed observation with the much 
larger area experiment on HEAO-1 A2 was made, but one burst was observed; the 
data on this one burst provides substantially more Information about the 
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spectral evolution than Is possible with the OSO-8 data. While the 
Infoimatlon available from OSO-8 on Individual bursts Is limited, these 
observations do give Information on the range of light curves and spectral 
properties. Comparison of results Indicates that the one burst observed with 
HEAO-1 A2 was similar to some of those observed with OSO-8. 

A. Rec urrence Time 

Figure 1 shows both the time that elapsed between bursts and the average 
steady flux from the source as a function of time for all the OSO-8 
observations. Gaps In the data of ~ 50% of an ~ 6000 s orbit mean that It Is 
quite common to miss bursts and that the burst Interval can appear to be 
multiples of the true value; from Figure 1 it Is relatively easy to see when 
this occurs. Removing those points that reflect missed bursts gives a mean 
recurrence Interval of ~ 16 ks (0.14d). In Figure 1 averages of are 
shown for several shorter Intervals Indicating drifts in the recurrence time 
from 13 ks to 18 ks. This is not unusual since jitters of up to ~ 10% in the 
recurrence times are found to be typical of most X-ray bursters (Lewin 
1977). 

The steady flux Is subject to major dips approximately every 50 
minutes. Sometimes less deep anomolous dips occur between major dips (Paper 
1; Walter et al. 1982). In Figure 2 we show some of the OSO-8 observations 
with a time resolution of 82s. It can be seen that in addition to the dips 
there are variations In the peak non-burst flux over timescales on the order 
of a day. These variations are also evident in the average count rates shown 
in Figure ?. The errors in Figure 1 reflect the dispersion In values for the 
non-dip phases. While the observed count rate varies by a factor of ~ 2, It 
is also associated with spectral changes (Paper 1) that give a much smaller 
change in the total 1-60 keV flux Fj^ from 1.0 to 1.25 x 10“® ergs cm"^ s"^. 
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(Hers and In the following the subscript “s" Is used to denote the "steady" 
source properties and "a" the apparent quantities measured by the observer.) 
There Is some evidence for an anti -correlation of this flux with the average 
Interval between bursts; when the flux Is low, the recurrence time Is 
longer. It Is not possible to make a meaningful test of more detailed 
correspondence with the available coverage, given the disruption by the dips. 

Although the dips were occasionally deep, they usually did not reduce the 
flux to zero and bursts were observed at the edges of dips. Figure 3 shows a 
histogram of the phases of the 50 minute period at which the bursts occurred 
where phase zero Is defined as the deepest dips. The distribution appears 
random, as would be expected If the occul tat Ions and the bursts were 
Independent. The data are not adequate to prove that the burst source Is 
occulted, although no bright burst peak occurred at the center of a deep dip. 
B. The HEAO-1 A2 Burst 

As reported In Becker et al. (1977), the bursts from MXBl 91 6-053 have 
e- folding decay times as short as 5 s. For the burst observed with HEAO-1 A2 
the count rates were high enough to derive useful measures of the source size 
and temperature on time scales short compared to this. In addition, the burst 
could be followed with less precision for an additional ~ 100 s. Figure 4 
(upper) shows the whole HEAO-1 A2 observation by the HED (2-60 keV) and Figure 
4 (lower) shows light curves In both the HED and the MED (2-20 keV) for the 
first 40 s of the burst. One short (30 s) low point In Figure 4 (upper) 
stands out as a candidate for a small absorption dip. The time 2 x 50 min. 
before that Is a few minutes after the burst and there Is no notable evidence 
of a dip. The deepest dip here Is very shallow and short in comparison to 
many observed In the OSO-8 data (see Fig. 2). Thus this burst was probably 
not affected by a dip. 
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Because of the Ifmited telemetry rate, spectral resolution was traded for 
temporal Information: at 80 ms, the shortest time scale available, we have 2 
broad band "colors'*. Additional bands are available every 1.28 s. We have 
used eight 1.28 s colors, 4 from the argon MED (Medium Energy Detector) and ^ 
from the xenon HED 3 (third High Energy Detector). Pulse height Information 
was available In 64 channels every 10.24 s. During the tall of the burst, 
when the changes are slow compared to 10 s and when the excess over the 
preburst spectrum Is predominantly of low energies, the MED pulse height data 
give the best determination of the energy spectrum. When the changes are fast 
compared to 10 s we can still use the 1.28 s colors to fit simple 3 parameter 
models. The consistency of the results was checked by comparing the spectral 
fits from the colors with the 10 s pulse height data. Here and In the 
following. In order to facilitate comparison with other burst sources, black 
body distributions were fitted to the burst spectra. For these fits, the 
column densities were about the same as required for the non-burst spectra. 

If we fix the column density, we can use the 80 ms data to determine a 
temperature and normalization. For the burst spectra In all cases we use the 
steady ^lux before the burst as background. Similar analysis procedures were 
used for a burst from MXB 1728-34 (Hoffman et al. 1979). 

The first part of the burst Is shown In Figure 5 with 80 ms resolution; 
the 2 colors have average energies of 5 and 11 keV, respectively. This time 
scale resolves the rise time to be 0.3 s with a possible lag of as much as 
100 ms of the high behind the low energy band. There Is a suggestion of 
temporal structure on time scales less than 1 s during the peak. Only during 
the first 4 s of the burst were the fits of the 1.28 s data to a black body 

model significantly better than to an absorbed thin thermal -bremsstrahlung. 

2 

Even for the black body model, the reduced x was 3 for the 1.28 s fits at 
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peak Intensity. When the Integrated Incident spectra Implied by the black 
body ^Its were tested as a model for the high energy resolution 10.24s data, 
the X values were acceptable at a level of confidence of 5% and 1%, 
respectively, for the two detectors with ^<5 channel pulse height 
Information. Estimates of the Integrated Incident spectra were generated and 
Figure 6 shows a comparison with the models as the burst evolvr^d. There are 
slight discrepancies between model and data, similar In both detectors, which 
suggest that the source Is Indeed not exactly black. But while the blackbody 
model Is not unique In fitting the data and the fits are not always 
acceptable, our data do not identify a systematic and reproducible 
discrepancy. 

The results of the black body fits are shown In Figure 7. The 
quantity e refers to a possible energy- Independent emissivlty such that the 
differential energy flux Is (R^/d^) e ttB where B Is a Planck spectrum, 

d, or djQ In units of 10 kpc. Is the distance to the source and ttR^ /d the 
subtended solid angle. The behavior of the temperature on timescales of 1.28s 
Indicated In Figure 7a Is typical of X-ray bursts. However, the Increase In 
the apparent radius by more than a factor of 3 over the 10 s decay of the 
burst Is unusuiil. Although Initial rises and declines are common (cf. Oda 
19811 the radius usually levels out during the cooling phase* In this burst 
It appears to gradually decline during the extended tall. Figure 7d,e, and f 
are based on the 80 ms data and are shown with various time resolutions. The 
(Independent) data for the 80 ms, 1.28 s and 10.24 s time scales give results 
which are In reasonable agreement. The rise 1s associated with a temperature 
Increase. While the uncertainties are large, ther§ Is a suggestion of 
temperature variations across the peak. During the first 3-4 seconds of the 
burst, no significant radius changes occur. The radius Increase starts when 
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the temperature starts to decline. In this respect these changes are 
typical; radius changes are usually associated with changes In temperature In 
the black bo4y model (Mhich would be expected If the model Is not correct as 
well as If the radius of the effective photosphere really changes). 

C. Burst Morphology 

Only two types of OSO-8 data were obtained: (1) 160 ms rates for the 
entire 2-60 keV band of the detector, obtained for four ~ 0.8 s Intervals 
during each satellite spin of ~ 10 s and (2) 64 channel pulse height 
Information accumulated over 41 s. The 160 ms data are shown for four 
"typical" bursts In Figure 8 (see also Becker et al. 1977 for e burst from the 
1976 observation.) While coverage of any one burst Is Incomplete, even among 
these chopped light curves there appear to be a variety of types among the 
1978 bursts: (a) bursts with peak fluxes of 0.7-2. 0 counts cm“^ s“^, rise 
times < 1 sec and decays which appear smooth with time constants of ~ 5s, (b) 
bursts similar to a^but lasting only 2-3s, (c) bursts with a plateau level 
of ~ 0.3-0. 7 counts cm“^ s"^ lasting 7-lOs with rise times which on occasion 
appeared as long as 2s, and (d) a combination of ^and £, that Is, a fast 
sharp peak and a slower decline from a plateau level. There does not appear 
to be a clear correspondence of these different burst types with steady flux, 
observed recurrence time, or phase with respect to absorption dips. The 
variety In the structure of the bursts has been noted In Hakucho observations 
of several other burst sources (Matsuoka 1981). While the behavior Is 
reminiscent of that of the transient 1608-522 (Murakami et al. 1980), the 
significant correlations between type and persistent flux or type and 
variation In peak luminosities were not present. 

The data averaged over a burst can be characterized by an average 
apparent kT^ that ranged from 1.5 keV (type £) to 3.4 keV (type_b), consistent 
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with there being a correlation between Intensity and temperature. However an 
average kT^ of 2.3 keV was obtafned for one type £ burst and SAS-3 results 
show a variety of peak temperatures for any particular peak Intensity range 
(Comlnsky 1981). If the response of the OSO“8 detector Is considered, the 
peak Intensities were between 1.7 and 4.8 x 10"® erg cm"^ s*^ (for peaks 
between 0.7 and 2 cts cm"^ s“^) and peak kT^ In the range 2. 6-3.0 keV. For 
the same ranges the Implied apparent radii of a spherical black bot(y were 
(6±2) X d^Q km, In the range of the radii seen In the HEAO-1 A2 
observations. 

For comparison, the burst as seen with the HEAO-1 A2 MED Is plotted In 
Figure 8e with a time resolution of 160 ms and the same scales as In Figure 8 
a-d. It seems to be Intermediate between types c and d, with a sharp rise and 
a high temperature peak, but a relatively slow decay. The peak rate of type b 
bursts seen with OSO-8 and an average temperature of 3.4 keV Imply an apparent 
radius of 4.5 x d^o similar to that at the beginning of the burst seen 
with HEAO-1 A2 (Figure 7e). The rate for the very low temperature type £ 
burst Implies IT x djo results suggest a recurring correspondence 

of high temperature peak with “small” radius and low temperature plateau with 
“large” radius. 

III. Derived Burst Quantities 
A. Ratio of Steady to Average Burst Luminosity 

The quantity o « f^sa^r^^^ba where F^^ and are the apparent 
steady and burst fluxes, reflects the nuclear energy liberated In the flash 
from the matter which has fallen onto a neutron star. Accounting for the 
observed ratios of ~ 100 was an outstanding success of the Helium flash model 
(Joss 1978; van Paradijs et al. 1979). This ratio does depend on parameters 
like the neutron star mass, the type of fuel. Its accretion rate and the 


conditions under which It accumulates (e.g. Team 1981a, b; Ay^sll and Joss 
1982), although It Is not yet a clear discriminant. For HXB1916-053 the OSO-8 
observations give F^^ 1.7 x 10’® ergs cm’^, approximately Independent of 

the changes In F,^ because of the anticorrelatlon with x^. For the burst seen 
by HEAO-1 A2 dt ■ 1.4 x 10’^ erg cm’^. This quantity also tended to stay 
constant, Independent of the burst types; I.e. there Is an Inverse correlation 
of peak flux and duration. Thus during both OSO-8 observations the values 
of a averaged over about one day appeared to stay constant at ~ 120. 

B. Peak l.um1no$1t1es 

HEAO-1 A2 observed a maximum of 1.8 x 10*® ergs cm“^s*^ averaged over 160 
ms. For OSO-8, peaks of > 4 x 10’® ergs cm’^s’^ are Indicated. The apparent 
luminosities, assuming a spherically symmetric source are (2-5) x 
10®® djQ ergs s’^ Implying luminosities at the surface of the neutron star of 
(6 - 14} s d|g/(3 g^) where S^o • 1 - 2 GM/Re^* (1 + Z)'* takes the 
mlivimum value ~ 1/3 for neutron star models (e.g. f: ■ 1.41 R « 6.5 km) 

(see Arnett and Bowers 1977). Considering that at the surface 

■•ED ' *00 TTfxr ‘nr’ * • 

where X Is the fractional abundances by mass of H, the ratio of the peak to 
Eddington limiting flux was 

^ - (1. 4-3.2) (3 (1+X)dj(, (1) 

For a 1.4 Mq neutron star at 10 kpc, accreting solar abundances, this takes 
values 2-4, and 1-2 for accretion of He. In the spherical blackbody model 
(neglecting effects of magnetic field or emissivlty appreciably less than 
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unity) the Instantaneous ratio of flux to Eddington limiting flux should be, 
for a neutron star with R > 3 6M/c^, 





^ <^9oo> 


(R/6.5 km)^ 
0 




(2) 


As lioldman (1979) pointed out, for kT^ ~ 2-3 keV this exceeds unity unless R 
and M are less than expected for neutron stars. We call a “standard'' neutron 
star one with mass and radius In the range discussed by Arnett and Bowers 
(1977) (In particular, with an equation of state harder than that used for 
curve B of their Figure 11). Marshall (1982) has pointed out that under these 
conditions the Eddington limit Is exceeded as long as kTg > 1.6$ keV for X = 1 
and > 1.96 keV for accreting He. Conversely If a standard neutron star were 
black, kTjj 3 keV would Imply > 9 times the Eddington limiting flux for solar 
abundances. We note however that this Interpretation would Imply for the 
HEAO-1 A2 burst d^Q » 2.3 (3goo)"^/^ (R/6.5 km), which Is always greater than 
2 for the allowed M a?id R. Thus for this source emissivlty or asymmetry Is 
Important, If the source Is not to be on the other side of our galaxy* 

Radius of a Spherical Black Body 
The apparent “radii", or more generally, e ^3/^10* 
data ranged from 4-15 km. Then for a sphere with R > 3 GM/c^, 


(3 9 


00 


)^^^(2,3 - 8.6) -iL km 


(3) 


The smaller value refers to the first 3s of the burst, the large value to the 
average high values 20-40 s after the onset. 
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ft. Accretion Rate 


The steady flux arises from mass accretion at a rate M at the neutron 

• 1/2 * 

Star surface which appears from afar to b » M ■ 9qo energy 

released Is 


Hc^ - 1] 


(4) 


with * cL appearing In the X-ray band. Here 

•-X • Soo'^ '^sa 


(5) 


where the correction factor f represents asymmetry In the emission of effects 
of limb darkening or obscuration {cf. Mllgrom 1978) which might be expected 
for the high Inclination Implied by the absorption dips. Thus the observed 
stea<(y flux Implies 


*^a " 4.7 x 10“^° ^ d^Q yr’^, 

^ “ ^00 


( 6 ) 


If the steacly emission Is not beamed and f < 1 


> 4.7 x 10“^° d^Q yr"^ 


(7) 


IV. EVOLUTIONARY CONSIDERATIONS 

In the following we assume that the mass losing star approximately fills 
Its Roche lobe. The combination of the Roche lobe relation for mass ratios 
less than 0.8 and Kepler's third law yield the well known relation P /p = 3.78 
X 10^ (Faulkner, Flannery and Warner 1972), where P and p represent the 
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orbital period and the mean density of the mass losing component 
respectively. For a 50 minute orbital period the secondary's mean density Is 
159 gm cm"^. In such a tight binary system a lower limit to the mass exchange 
rate can be estimated from the loss of angular momentum via gravitational 
radiation (GR) (Faulkner 1971; Paczynski 1967). If the evolution driven by 6R 
(or evolution driven by some other mechanism) Is rapid compared to the thermal 
time scale of the mass losing star, thermal disequilibrium effects will be 
Important (Paczynski 1981). It Is convenient to separate the dlscusslc*^ Into 
the cases where such thermal effects are not (case I) or are (case II) 
Important. The results are summarized In Table 1. 

Case (I) -- Thermal Equilibrium 

If thermal disequilibrium effects are not Important, one can adopt 
equillbrum mass-radius relations. As stated In Paper I, possible solutions 
Include (1) a 1.11 He main sequence star (based upon models calculated by 
Paczynski 1971), (2) a 0.044 “hydrogen" (X = 0.75) white dwarf and (3) a 
0.009 Mq He white dwarf (the latter two solutions based on the models 
calculated by Vila 1971). No self consistent solution for a hydrogen main 
sequence star was given In Paper I, since the best estimates of the mass 
radius relation for the lower main sequence Imply a mass (0.076 M^) below the 
minimum for hydrogen burning (0.085 M 0 according to Graboske and Grossman 
1971). The difference Is small however and subtleties In the treatment of the 
equation of state might affect the calculated endpoint of the main sequence or 
the rfiass radius relation, so that a hydrogen burning companion may still be 
possible. Recently, D'Antona and MazzItelH (1982) have constructed low mass 
models assuming constant mass loss rates. 

For M £ 10" yr“ , (corresponding to models which are In thermal 
equilibrium) their models give the Inferred mean density at 0.093 and 0.10 



for Z - 0.02 and 0.001, respectively (where Z Is the metal abundance by 
mass) . 
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By generalizing the expression for mass transfer rates Induced by 
gravitational radiation losses, assuming conservative mass transfer (cf. 
Faulkner 1971) one obtains, for a 1.4 neutron star, mass transfer rates of 
6 X 10"®, 7 X 10"^^, and 3 x 10"^^ yr"^ for solutions (1), (2) and (3), 
respectively. For hydrogen burning main sequence mass-radius relation quoted 
by Whyte and Eggleton (1981) (extrapolated below 0.085 Mq) we find a mass 
exchange rate of 1 x 10"^® M 3 yr"^. If Instead we take R ~ M, (used by LI et 
a1. 1980), the mass exchange rate Is 80% higher. The models by D’Antona and 
MazzItelH for constant, low mass loss rates Imply ~ 2 x 10"^® M 3 yr“^ For 
the solutions with secondary masses < 0.1 M 3 , the mass transfer rate Is 
proportional to M^2/3^ where M^ Is the neutron star mass, so that a higher M^ 
could Increase the rate only about 50%. 

B* Case II — Thermal Disequilibrium 

Whenever the mass transfer proceeds on a time scale which Is short 
compared to the thermal time scale of the mass losing component, the thermal 
history of the latter will affect the stellar radius. In fact, a sizable 
departure from thermal equilibrium mass radius relations Is expected. RJW 
have recently Investigated the evolution of binary systems composed of main 
sequence stars and compact objects with gravitational radiation losses 
Included. As pointed out In RJW, the effect of thermal disequilibrium becomes 
pronounced as the system approaches Its minimum orbital period. A minimum 
period Is expected due to the transition In stellar structure from a hydrogen 
burning star to a degenerate one. If we adopt their case 6 which corresponds 
to a hydrogen depleted composition with X = 0.2 and Interpolate In their 
Figure 8 , we obtain a mass of about 0.04 M 3 for the mass losing star. The 
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corresponding mass transfer rate near the minimum period of 46.5 minutes Is 
approximately 4 x 10"^^ yr"^. 

Departure from thermal equilibrium can also be Important for secondaries 
composed of a helium core and a hydrogen rich envelope. Such composite models 
were suggested by Faulkner, Flannery, and Warner (1972) as possible 
progenitors to HZ 29 AM CVn. For such evolved secondaries, the stellar 
envelope would contract once Its mass reached a critical minimum value (cf. 
Refsdal and Welgert 1971). Consequently, mass transfer would cease as the 
star no longer would fill Its Roche lobe. However, since the rate of stellar 
contraction decreases as the secondary evolves to the white dwarf stage, the 
loss of angular momentum by gravitational radiation would eventually reduce 
the size of the Roche lobe until mass transfer Is reinitiated. In such an 
evolutionary scheme a mass transferring system at the 50 min. orbital per1o-l 
Is a distinct possibility. We have computed such an evolutionary sequence for 
low mass composite secondaries (M « 0.11 and 0.15 M^) In order to test this 
conjecture. 

The Initial models are characterized by a Population I composition (X * 
0.70, Y » 0.28, and Z » 0.02) and were constructed and evolved utilizing the 
stellar structure program described In Eggleton (1971,1972). The assumptions 
and constitutive relations are described In Taam (1983). The two model stars 
were evolved at constant mass to the envelope contraction phase and then 
placed In a binary. The period and total binary mass corresponded to 50 min 
and 1.51 Mg and 58 min. and 1.55 Mg for sequences 1 and 2 respectively. The 
results for both sequences are summarized In Table 1. We find that mass 
transfer rates as high as 2 x 10"^® Mg yr"^ can be achieved at an orbital 
period of 48 min for sequence 1, while for sequence 2 the mass transfer rate 
Is lower (0.4 X 10“^® Mg yr“^). The mass transfer rates for each sequence are 
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plotted as a function of orbital period In Figure 9. Secondaries as massive 
as 0.15 Mg contract rapidly at these periods Implying low mass transfer 
rates. A third sequence (characterized by a secondary of 0.11 Mg and an 
Initial period of 71 min) was calculated to determine the effect of even lower 
mass secondaries at SO min. It Is found that the mass transfer promoted by 
gravitational radiation Is slower than for sequence 1 because the angular 
momentum losses are less efficient for smaller mass stars at a given orbital 
period (M ~ 0.052 Mg at a 48 min orbit). 

V. COMPARISON WITH THE OBSERVATIONS 

Assuming GR to be the dominant driver of mass exchange, one can constrain 
the distance for each of the models discussed In the previous section. Since 
the Inferred mass transfer rate Is > 4.7 x 10“^® d^lQ Mq yr’^ (see eq. 7), the 
source would have to be closer than 5 kpc for the solutions with M^ ^ lO’^^ Mg 
yr"^. For the He degenerate solution the limiting distance would be ~ 1 
kpc. For the H depleted (X « 0.2) solution In which thermal disequilibrium Is 
taken Into account in the evolution the distance is ~ 3 kpc. For the 
composite model GR would saturate the mass exchange at 6. 0-6. 5 kpc. 

The He main sequence solution was ruled out in paper I and by Walter et 
al. (1982) on the basis of the optical appearance of the source. The bursts 
provide another argument against it, If we assume the type I bursts are 
thermonucl ear runaways in freshly accreted layers on a neutron star (Joss 
1978; Taam and Picklum 1979). The solution Implies gravitational radiation 
driving the mass exchange at such a high rate (in excess of the rate for which 
the steady luminosity equals the Eddington luminosity) that the luminosity 
from the nuclear outburst would not be sufficient to dominate that from 
accretion. The He degenerate secondary solution (see above) Is also 
inconsistent with the burst data. In particular, the expected burst 
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recurrence tiroes would be roore than an order of magnitude longer than those 
observed. Thus all solutions involving a He secondary are very unlikely. 

Thermonuclear flash theory (Joss 1978; Taam 1981a, b, 1982; Ayasli and 
Joss 1982) predicts that the peak burst luminosity does not exceed the 
Eddington limit. Based upon equation (1) the peak luminosities of most bursts 
observed from the source would be equal to the Eddington limit at ~ 6 kpc. 
Thus, the composite model provides a self-consistent scenario if the bursts 
peak at about the Eddington limit. However, there is evidence that points to 
burst luminosities exceeding the Eddington limit by factors of 2-10 even if 
black body emission is not assumed (Inoue et al. 1980; van Paradijs 1981). 
Thus, comparison with other bursters would favor a distance near 10 kpc if the 
galactic center is 9 kpc away and not as close as 7 kpc (cf. Frenk and White 
1982; Glass and Feast 1982). The distribution of known burst sources about 
the galactic center is uno^rtain, but as Inoue et al. (1980) discussed, the 
group is unlikely to be clustered far from the galactic center. Thus the 
composite model solution is only marginally consistent with both 6R and 
results on other bursters. Furthermore we must bear in mind that to get this 
consistency we are forced to use the lower limit on the mass transfer rate. 

If the assumption is made of spherical black bodies- comparison of the 
luminosities of the bursters is nearly the same as comparison of the radii, 
since the temperatures of the i')lack body fits are typical. The values 
of ^ (4-15) dio km observed for bursts from MXBl 916-053 are near the 

range 5-12 km observed for other bursters (Inoue et al. 1980; van Paradijs 
1978) and would still overlap this range if djo “ 0*®* 

Ayasli and Joss (1982) and Taam ( 1 981 ab, 1982) find that for short 
recurrence times of less than 5 hours the hydrogen rich layers participate in 
the flash (see also Fujimoto, Hanawa and Hiyayi 1981) and the mass accretion 
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rates should exceed 5 x 10"^® yr"^, (although effects due to the evolving 

thermal state and compositional history of the neutron star envelope could 
reduce this rate somewhat). This Is In agreement with the observed steady 
flux If the source Is near 10 kpc. At this distance the bursts, then, are 
more typical In peak luminosity and Inferred radii to other burst sources. In 
addition nuclear flashes of this type (H-He) are consistent with a ~ 120 found 
for this source. 

VI. DISCUSSION 

If the burst emission Is spherically symmetric from the entire surface, 
flash theory (at least In the limit cycle approximation) Implies a mass 
transfer rate higher than 6R can drive from any of the types of secondary 
considered In SIV, except possibly the composite model (with M 0.1 M^j). A 
spherically symmetric stea((y source Is then at such a distance (10 kpc) that 
the burst luminosities are typical. 

For 4U1626-67, Joss, Avni and Rappaport (1978) and LI et al. (1980) 
concluded that GR was marginally sufficient to drive the mass exchange for a 
0.1 Mg hydrogen burning secondary because the 10"® Mg y"^ estimated from the 
spin up rate was very uncertain. For MXB1916-053/4U1915-05 the bursts 
constrain the exchange rate. In both cases the discrepancy with GR Is least 
for a secondary at the end of the main sequence and Is minimally a factor 
of ~ 2. This discrepancy Is not so much greater than the uncertainties In the 
arguments as to be unquestionable but a second source now shows the same 
discrepancy on the basis of quite different arguments. 

Several mechanisms have been appealed to for driving higher rates of mass 
transfer than seems possible with GR alone. Nuclear evolution (Faulkner 1974; 
Joss, Avni and Rappaport 1978; RJW) may be a likely possibility for more 
massive systems with periods longer than ~ 4*^. For lower mass systems the 
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mechanisms of (1) ^‘magnetic braking'* and (2) transport of angular momentum 
out of the systr;m with ejection of some material In the mass exchange process 
may be Important. 

Coronal magnetic field structures may force matter to corotate to large 
distances, so that If It then escapes, It carries off substantial angular 
momentum. Even If matter does not escape, the Interaction of the fields and 
corona with surrounding material may exert a torque. The angular momentum 
that could be carried away by modest mass loss (10“^^ Mq y‘M has long been 
recognized to be significant and can Influence binary evolution (Mochnacki 
1981 and references therein). Verbunt and Zwaan (1981) estimated that If the 
empirical .braking law of single stars can be carried over to the close binary 
situation, this should only dominate over GR for secondary masses above 0.3 
Mq. The limit can be reduced to about 0.1 (Taam 1983) when the effects due 
to thermal disequilibrium and the change In the moment of Inertia of the mass 
losing star are Included. Rappaport, Verbunt and Joss (1983) found that with 
different braking laws, magnetic braking would dominate for lower masses, but 
the minimum periods of those models were too long. While It seems probable 
that these calculations cannot adequately represent the binary situation, the 
effect seems likely to be Important for low mass non-degenerate and possibly 
magnetic degenerate companions. 

A large amount of angular momentum must be removed from the exchanged 
material before It accrete onto the neutron star and It seems plausible that 
some. If not a Targe fraction. Is ejected from the system. Paczynski (1967) 
considered an exchange of angular momentum with an outer ring of the disk. The 
angular momentum may return to the orbit (Papal oizou and Pringle 1977; Weber 
1979), but perhaps this could be prevented by effects not yet fully 
considered. Self-excited winds as a way of ejecting angular momentum was 
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noted by many authors and indeed Begelman, McKee and Shields (1983) have found 
that a substantial Mind is expected. The effect does not turn out to be 
substantial on the evolution of binaries Mith non-degenerate companions (Whyte 
and Eggleton 1980; Kieboom and Verbunt 1981; RJW), and ROW showed this should 
be true for any model in which the mass loss is proportional to the mass 
exchanged. However, as we have shown, for this source a rather modest factor 
of 2 may be all that Is needed. This mechanism might be even more Important 
for a degenerate mass losing star. Kieboom and Verbunt (1981) noted that for 
the case of standard abundances, if a small amount of matter carried off the 
angular momentum, an order of magnitude effect could result; the mass exchange 
for a 0.05 degenerate companion was as high as 10"® Mg yr"^ for ~ 10® yr. 

VII. SUMMARY 

We have used the burst properties of the 50 min binary MXBl 91 6-035 to 
estimate the mass transfer rate between the binary components and hence 
investigate the consistency of the GR evolutionary scenarios developed to 
date. Both degenerate and non-degenerate helium secondaries are incompatable 
with the burst properties and can be ruled out. If it is assumed that the 
peak burst flux is at least at the Eddington limit then a hydrogen degenerate 
companion will give a mass transfer rate at least a factor of 3 below that 
implied by the steady flux. On the other hand, the composite. Inhomogeneous 
model is marginally consistent with the GR driven mass transfer rate. Current 
flash theory prefers a rate larger by about a factor of 2 as does comparison 
of its burst properties to other sources. In that case there may be another, 
dominant mechanism for losing angular momentum from the system. Recent work 
Implies that If the angular momentum extracted from the accretion disk is lost 
from the system rather than fed back into the orbit, the transfer driven by GR 
could be sufficiently enhanced. Magnetic braking Is another possible 



Influence, although probably more Important for more massive secondaries. 
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POSSIBLE SOLUTIONS FOR THE SECONDARY 
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TYPE 

M(M^) 

Equilibrium 


H burning 

0.076 

Z « 0.02 

0.094 


0.098 

Z « 0.001 

0.103 

He burning 

1.111 

H degenerate (X * 

.75) 0.044 

He degenerate 

0.009 

Non-equi librium 


H burning (X » .2) 

-0.04 

Composite 

0.11 


0.15 

OBSERVED 



* Whyte and Eggleton (198C) 

^ LI et al. (1980) 

^ D’Antona and MazzItelH (1982) 

^ Paczynski (1971b) 

® Vila (1971) 

^ Rappaport, Joss and Webbink (1982) 
9 this paper 


R(Re) 

”gr<”o 

REFERENCE 

0.088 

9.78 X 10-11 

a 

0.094 

1.78 X 10*10 

b 

0.096 

1.70 X 10-10 

c 

0.097 

2.03 X 10*10 

c 

0.215 

6.1.) X lO’O 

d 

0.0/3 

6.96 X 10*11 

e 

0.043 

2.82 X 10*12 

e 

0.07 

~ 4.0 X 10*11 

f 

0.10 

1.8 X 10*10 

9 

0.11 

4.0 X 10*11 


— 

>4.7 X 10*10 j2^ 
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FIGURE CAPTIONS 


Figure 1 (upper) - Time since the last burst observed during the OSO-8 
observations (•). The lines are averages of Intervals between successive 
bursts. See text for definition of burst types a-d. (bottom) - Average 
non-dip count rate as a function of time. 

Figure 2 - Counts cm“^s“^ averaged over 82s for 2 Intervals 1 d long of 
OSO-8 observations from 1978 March. Typically, every other dip was obscured 
by earth occultatlon. The average non-dip count rates were variable by almost 
a factor of 2. Altogether ~ 50 dips were observed. 

Figure 3 - Histogram of number of bursts falling In phase Intervals of 0.1 
relative to the phase of dip center for 1976 (hatched) and for 1978 (open). 

Figure 4 (upper) - Light curve of 4U1 916-05 with the HEAO-1 A2 HED, similar 
to the detector used In the OSO-8 observations but with 2.9 times the area. 
Arrows point out a dip and the times when the same phase of a 50 min period 
was covered, (lower) - The HED and MED burst light curves with 1.28 s 
resolution for the first 40s of the burst. 

Figure 5 - The first 8 s of a burst are shown with 80 ms resolution In the 
MED (argon) detector (upper) for which the average energy Is 5 kev and the 
second layer of the HE03 (xenon) detector (lower) with average energy of 11 
keV. 
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Figure 6 - Average incident spectra for three Intervals of the burst seen by 
HEAO-1 A2. 

Figure 7 - The spectral parameters derived for the HEAO-1 A2 burst, (a) kT^ 
obtained from 1.28 s color fits (solid crosses) and 10.24 s pulse height data 
(dashed crosses) (b) V^IO <c) flux (1-60 keV) for the 

unabsorbed black botly spectral model; the solid line Is the Integrated energy 
per unit area received from the burst, for which the scale Is 10“^ ergs 
cm‘^. (d) kT^ for the beginning of the burst from the 80 ms data: 320 ms 

averages (dashed crosses) and shorter averages (solid crosses) along with the 

1.28 s results ( ). (e) ^a/^10 same as (d). (f) Flux for 

the same as (d). 

Figure 8 - Light curves of 50 s of 4 OSO-8 bursts and the 1 HEAO-1 A2 burst 
with 160 ms time resolution, (a), (b), (c) and (d) refer to the 4 types 
discussed In the text. The HEAO-1 A2 light curve Is from the MED. 

Figure 9 - Mass exchange rate driven by GR as a function of the binary 
period for three sequences. The secondary mass (In Mg) is Indicated at 
selected times during the evolution. 
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